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Dark matter particles can be captured by the sun with rates that depend on the dark matter 
mass and the DM-nucleon cross section. However, for masses below ~ 3.3 GeV, the captured dark 
matter particles evaporate, leading to an equilibrium where the rate of captured particles is equal 
to the rate of evaporating ones. Unlike dark matter particles from the halo, the evaporating dark 
matter particles have velocities that are not limited to values below the escape velocity of the galaxy. 
Despite the fact that high velocities are exponentially suppressed, I demonstrate here that current 
underground detectors have the possibility to probe/constrain low dark matter parameter space by 
(not)-observing the high energy tail of the evaporating dark matter particles from the sun. I also 
show that the functional form of the differential rate of counts with respect to the recoil energy 


in earth based detectors can identify precisely the 
particle in this case. 

Despite the overwhelming amount of evidence from 
galaxy structure formation and cosmology in favor of the 
existence of dark matter (DM), so far experiments have 
failed to conclusively detect directly or indirectly DM. 

A lot of scientific effort and resources have been allo¬ 
cated in order to detect DM. In particular, underground 
detectors have imposed strict limits on the type of DM 
particles that can be viable, constraining the parameter 
space of the mass and the cross section of DM interacting 
with baryons. The basic principle of direct detection is 
simple. A DM particle interacts with a nucleus of the de¬ 
tector deposing an amount of energy that is detectable. 
Different experiments apply different techniques on how 
they observe the recoil. However, direct detection rates 
have limitations. Obviously a small DM-nucleus cross 
section reduces the probability of interaction. Similarly, 
heavy DM particles have low number densities and con¬ 
sequently low flux. Low DM masses are also difficult 
to probe simply because the DM particle does not have 
enough energy to trigger the detector. This is true re¬ 
gardless the exposure that an experiment can achieve. 

No matter what is the velocity distribution of DM parti¬ 
cles in the galaxy, their velocities are below the escape ve¬ 
locity of the galaxy. Therefore below a sufficiently small 
DM mass and a given detector energy threshold, no DM 
particle can be detected. This is the reason why low DM 
masses are not probed by direct DM searches. 

However, as I will demonstrate in this paper, it is pos¬ 
sible to probe lighter DM masses with current detectors 
and energy thresholds due to a flux of DM particles that 
after being captured by the sun, leak out via evaporation. 
These particles can arrive in earth with energies that are 
high enough to produce a detectable recoil. Therefore 
not only it is possible to probe relatively lighter DM, but 
additionally the spectrum has features that distinguish 
it clearly from heavier DM candidates. The possibility of 
detecting DM particles that have been captured by the 
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mass and the cross section of the dark matter 


sun, on the earth, has been studied in the past [1, 2]. In 
these two seminal papers, the flux of DM particles bound 
to the sun and having orbits that can reach the earth was 
estimated. Damour and Krauss considered particles that 
have been captured by the outer layers of the sun and due 
to perturbations from other planets, the orbits evolved 
to eliptical ones that do not cross the sun anymore and 
therefore do not lose further energy. Particles in these or¬ 
bits can accumulate for billions of years and since some 
of the orbits can reach the earth, these particles are po¬ 
tentially detectable. Additionally, this scenario has been 
studied numerically [3, 4]. It should be emphasized here 
that this present paper studies a fundamentally different 
scenario. Instead of looking at loosely bound DM parti¬ 
cles that have orbits that cross the earth, I focus on light 
particles that have had the time to thermalize with nu¬ 
clear matter inside the sun. The tail of the distribution 
of these particles corresponds to velocities above the es¬ 
cape velocity of the sun and this is exactly the spectrum 
of particles that I consider here. 

Generally, the number of DM particles in the sun N is 
determined by 

^=F-CeN-CaN^, ( 1 ) 

at 

where F is the capture rate, and Ce,a are coefficients 
related to the evaporation and annihilation of DM re¬ 
spectively. From the above equation it is clear that if 
Cg >> CaF, and C~^ is much smaller than the age of 
the solar system, evaporation dominates the whole pro¬ 
cess and an effective equilibrium between the accretion 
rate F and evaporation has been established by now. In 
other words, if the above conditions are satisfied, DM 
particles leak out from the sun with the same rate as they 
are captured. Let us first estimate the capture rate [5-8] . 
The capture rate is 

where Mq and Rq are the mass and the radius of the sun, 
vq the velocity dispersion of DM in our galaxy, pum and 
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the local DM density and the DM mass respectively. 
The sum runs over all different chemical elements present 
in the sun. I am going to consider for simplicity only hy¬ 
drogen and helium. Ei is the maximum energy per DM 
mass that can lead to a capture due to a collision of DM 
with element i and is given by Ei = 'yiGMQ/R q (1 — 7^) 
where 7^ = -I- mMiY is the average frac¬ 

tion of energy that the DM particle loses after collid¬ 
ing with a nucleus Ni. Note that for energy larger than 
Ei, even if the particle scatters inside the sun will lose 
on average an amount of energy that is not sufficient 
to bind the particle gravitationally to the sun. Finally 
fi represents the probability that scattering will take 
place, fp — 0.896pfjp/fJcrii5 if ^p^pj^crit ^ 1 and 1 
if i-pOpjacrit > 1; where is the mass fraction of hy¬ 
drogen in the sun that is taken to be 0.75. For he¬ 
lium, jue = 0.89 X 4e_f/eO-p/crcr*i, if ^^Hef^pIcTcrit < 1 
and 1 if ^euecrplUcnt > 1, where tHe = 0.24. acnt = 
mpR^/MQ ~ 4 X 10“^®cm^ is roughly speaking the cross 
section above which every particle that will cross the 
sun will scatter. Note that everything is expressed in 
terms of the DM-proton cross section ap. I consider spin- 
independent interactions and therefore the DM-helium 
cross section will be ane = where cor¬ 

responds to the reduced mass of DM with nucleus i and 
xl = 4 for helium. Since I am interested in low DM mass, 
tJHe — 16(Jp. 

Now let us focus on evaporation. This effect has been 
studied extensively in the case of the sun [9-13]. Unless 
one assumes unreasonably high DM annihilation cross 
section, it has been shown [12, 13] that for DM masses 
below ~ 3.3 GeV, DM particles get effectively evaporated 
out of the sun. In this case the steady state solution of 
Eq. (1) will give an equilibrium between captured and 
evaporated DM particles. Although the exact formula of 
Ce has been estimated [12, 13], it will not be needed here. 
As long as I consider particles below 3.3 GeV, the rate of 
evaporation will be equal to that of capture. Therefore 
the overall evaporation rate will be given by Eq. (2). Let 
us now determine the spectrum of the evaporating DM 
particles. In general there are two possibilities. If the 
DM-nucleon cross section is large and the mean free path 
of the DM particle small, the captured population of DM 
will thermalize fast with nuclei and the DM distribution 
will be a Maxwell-Boltzmann one with a DM temperature 
equal to the one of the star at a particular position. How¬ 
ever, if the DM-nucleon cross section is small, captured 
DM interacts over several orbits and therefore there is no 
single temperature that picks up. The distribution is not 
a Maxwell-Boltzmann one, but it can be approximated 
as one although the DM “effective” temperature is differ¬ 
ent from that of the star [9, 13]. This approximate dis¬ 
tribution should look like f{v,r) ^ exp[—E{v,r)/kT^], 
where is the “effective” temperature of the DM, and 
E{v,r) = m^v‘^12 + m^V{r) is the total energy, V{r) 
being the gravitational potential as a function of r inside 
the sun. I would like to estimate the spectrum of evap¬ 
orating DM particles now. From this point of view, the 


details of the spatial dependence of the distribution are 
irrelevant since I consider particles that are at distance 
Rq from the center of the sun with a velocity higher than 
the escape velocity of the sun. Therefore the spectrum 
of evaporating DM particles will be given by 

m 

f{v) = Ae , V > Vs, (3) 


where A is a constant to be determined and Vs = 
{2 GMq/RqY/"^ is the escape velocity from the surface 
of the sun. However, this spectrum of evaporating DM 
particles does not remain the same when the DM parti¬ 
cles arrive on earth. Let us find the spectrum of velocities 
at the earth f{r,v) by using the collisionless Boltzmann 
equation 


dt dr dvr 


(4) 


Due to the isotropy of the problem, df jdO = df /dcj) = 0. 
I am interested in a steady state solution and therefore 
df /dt = 0. Fr/rrix = —GMq/t'^ + GM^/{i — r)^ is the 
force due to gravity from the sun and the earth, where 
r is the distance from the center of the sun, is the 
mass of the earth, and i is the distance between the sun 
and the earth. Note that Eg = = 0. The generic 

solution of Eq. (4) is f{v/ — 2GMq/t — 2GMi^/[l — r)). 
This solution should match the boundary distribution at 
the surface of the sun f{RQ,v) given by Eq. (3). Upon 
using this boundary condition, the distribution in earth 
is 

f{i,v)=Ae (*' V > Ve, (5) 

where Ve = (flGMQ/l -|- 2 GM^/Ri^Y/"^ and I?® is the 
radius of the earth. I have omitted negligible terms 
of the order 0{GM(^/(). Note here that in Eq. (5), 
v'^ = v/. + where the values of vg and are 

equal to the ones at the boundary of r = Rq of Eq. (3) , 
since Eq. (4) does not involve derivatives of them. DM 
particles that evaporated from the sun arrive in the earth 
almost radially. This is because the ratio V 0 /vr of a par¬ 
ticle arriving in the earth varies from 0 to a maximum 
value of Rq/(- — 0.0046. The total flux of evaporating 
DM particles arriving in the earth is 

nOO p 


Recall that dfv = v^dvd cos Odcj) and as mentioned above 
the solid angle integral part does not extend to the full 
47r but it is constrained to the value mentioned above. 
The product of A with the angular integration part is 
determined by Eq. (6), thus leading to the following flux 
of evaporating DM on earth 


= Ge~ 


(7) 


where the constant G is 


G = 


F 



/ 2 I 2 2 \ 
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•^dv 


-1 


( 8 ) 
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Let us now estimate the number of counts registered 
in an underground detector taking into account both the 


flux of evaporating DM and regular DM halo particles. 
The differential rate of counts per recoil energy is 


dR 

dEft 


N'j' 



da 

dEn 


Ei{v)dv 


Px 



f{v)vd^v , 


(9) 


where Nt is the number of targets in the detector, 
= 0.3GeV/cm^ is the local DM density, Vesc = 
550km/sec is the escape velocity of our galaxy and 
Vmin = x/'> tT'nEh/2p? is the minimum velocity required 
to produce a recoil E^ in a DM collision with a nu¬ 
cleus of mass niN [p being the reduced mass between 
DM and nucleus). For the distribution f{v), a trun¬ 
cated Maxwell-Boltzmann function up to Vesc of the form 
f{v) = A/’exp[—(u -|- Vh)'^/ vq\ is used, where A/" is a nor¬ 
malization constant [14, 15]. Vb is the velocity of the 
earth with respect to the rest frame of the halo. The 
value used here is Vb = (232 -|- 0.489 • 30) km/sec, which 
is the velocity of the solar system plus the rotational ve¬ 
locity of the earth around the sun (when the latter aligns 
maximally with the former). This value represents the 
best possible scenario for detecting halo DM particles. 
Note that the first term in the right hand side of Eq.(9) 
corresponds to the evaporating DM particles that arrive 
in earth from the sun, while the second one is the usual 
rate from incoming DM halo particles. A crucial obser¬ 
vation here is that although DM halo particles have an 
upper velocity of Vesc, the evaporating ones can have any 


energy by paying a price in an exponential suppression 
in the density. However this fact has important conse¬ 
quences because for a given threshold in DM detectors, 
and since v < Vesc, there is a mass below which DM 
particles from the halo can never have energies that can 
trigger the detector no matter how large the exposure is. 
On the contrary, for DM particles that have been cap¬ 
tured first by the sun and later evaporated, it is probable 
to detect the tail of their distribution since there is no 
upper velocity, if enough exposure is achieved. 

So far the value of the “effective” DM temperature 
has not been specified. This has been estimated in 
e.g. [9, 13] and it depends on the DM mass. However 
both papers gave results for DM masses above ~ 2 GeV. 
In order to find in much smaller masses of interest, 
I implement the method presented in [9]. Although as 
mentioned earlier, the actual distribution of captured 
DM is not an exact Maxwell-Boltzmann, it can be ap¬ 
proximated by such with a temperature T^. can be 
estimated by demanding no net flow of energy from the 
nuclei of the sun to the DM particles once a steady state 
has been achieved. The condition can be written [9] as 


d^r np{r) j d^uexp ( 


d^ 


Vp exp 


-TUpV^ 

2kT(r) 


ap\v^ - FpKAA) = 0, 


( 10 ) 


where np{r) and T{r) are the number density of nuclei 
and the temperature of the star at radius r respectively, 
and Vp are the velocities of DM and nuclei, and E is 
the total energy of DM (i.e. kinetic plus potential). Cp 
is the DM-proton cross section and (AA) is the energy 
exchange in a DM-nucleon collision. For simplicity I use 
a polytropic model of n = 3 as an approximation for 


the sun. If (j>{^) is the solution of the n = 3 Lane-Emden 
equation, np(^) = np{0)(j}{^)^, T(^) = Tc(j){^) and H(^) = 
{4:kTc/mp)[l — <(>(0]! where ^ represents a dimensionless 
radius defined as ^ Ci = 6.8968486 being 

the first zero of Eq. (10) can be rewritten in terms 
of the dimensionless quantities t = T^/Tc {Tc being the 
core temperature of the sun) and u = m^/nip as 



exp 


Ah' 

T 


T + 


[t - <?i(^)]C^d^ = 0 . 


( 11 ) 


Apart from a minor typo, the above equation is the same 


as that derived in [9]. For every considered DM mass. 
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I have solved numerically Eq. (11) in order to find the 
corresponding T^. 

cpd/(Kg.01keVee) 



FIG. 1: Number of counts per 0.1 keVee recoil energy nor¬ 
malized to an exposure of 1 Kg-day for a Si detector such 
the one used by DAMIC for DM-proton cross section of 
Gp = 10“®® cm®. A flat efficiency of 0.17 for the detector 
at low masses has been used. The dashed line corresponds 
to halo DM counts of DM with a mass of 1 GeV. The solid 
lines correspond to counts coming from DM evaporating from 
the sun with masses 1 GeV (thick solid), 100 MeV (medium 
solid) and 10 MeV (thin solid). Note that for DM mass of 10 
and 100 MeV, halo DM particles have not sufficient energy to 
recoil within the range shown in the plot. The vertical line 
represents the current threshold of DAMIG i.e. 40 eVee. 

Fig. 1 shows some characteristic results that demon¬ 
strate why the spectrum of evaporating DM particles can 
enhance the chances for direct detection of light DM par¬ 
ticles. It shows the rate of counts per recoil energy in bins 
of 0.1 keV, normalized to an exposure of 1 Kg-day, for a 
Si detector like the one used in DAMIC. I quote DAMIC 
here because it is one of the experiments with the lowest 
recoil energy threshold. I have used a flat efficiency of 
0.17 for the detector, deduced from [16]. Fig. 1 corre¬ 
sponds to DM-proton cross section of cTp = 10“®® cm®. 
For smaller cross sections, the rate of counts can be ob¬ 
tained by scaling the evaporating lines by (cTp/use)^ and 
the halo one by tTp/crae, where = 10“®® cm®. For the 
evaporating particles, one power of dp comes from the 
capture rate in the sun and one from the detection in the 
earth. I assume contact spin-independent interactions 
(using a Helm form factor) in this paper, although the 
generalization for spin-dependent is trivial. In the plot 
three cases of DM masses are shown: 1 GeV, 100 MeV 
and 10 MeV. One can see that halo DM particles roughly 
below 1 GeV cannot be detected. The reason is that light 
DM particles do not have enough energy to create recoil 
energies above the threshold. Note that the number of 
counts for the halo DM in the plot, correspond to the 
best detection scenario, i.e. the period of the year where 
the earth rotational velocity around the sun aligns maxi¬ 
mally with the sun velocity in the rest frame of the halo. 
On the contrary, light DM that has been captured by the 


sun and has evaporated after thermalization, have no up¬ 
per bound (apart from the speed of light). Although high 
velocities are exponentially suppressed, depending on the 
detector exposure, such particles can be detected. Fig. 
1 shows that slightly lower thresholds in direct detection 
can set limits down to DM masses of 10 MeV. This is 
impossible for halo DM particles that have no chance to 
be detected even with a threshold of order of eV. 

The spectrum of evaporating DM particles not only 
can probe/constrain parameter space that is not acces¬ 
sible by observing halo DM particles, but it can po¬ 
tentially determine accurately the DM mass and DM- 
nucleon cross section. As it can be seen in Fig. 1 for 
the case of m = 1 GeV, there is a change in the shape 
of the rate of counts in the detector. For low recoil en¬ 
ergies, halo DM particles dominate the counts. However 
the counts produced by halo DM particles drop sharply 
at the maximum recoil energy 'ysi’^^'i'^esc+Vb)'^ (jsi refers 
to the nucleus of Si and it is dehned below Eq. 2). Above 
this specific value of recoil energy only evaporating DM 
particles contribute to the counts. Therefore this drop in 
the number of counts per energy can lead to the exact 
determination of the DM mass and the DM-nucleon cross 
section. 

Some comments are in order here. One should make 
sure that there is enough time for the DM particles to 
thermalize with the interior of the sun. The issue has 
been addressed in [7] where it was shown that the char¬ 
acteristic time scales are of order of a year or smaller (for 
the bulk of the DM orbits). The second comment is re¬ 
lated to DM annihilation. The spectrum of evaporating 
DM particles shown here is valid whether one considers 
asymmetric DM or thermally produced symmetric DM. 
As it was argued in [13], for a DM mass of 3 GeV, the 
evaporation rate is equal to the annihilation one (with an 
annihilation cross section that of the weak interactions). 
For every 0.3 GeV below that mass value, the annihila¬ 
tion signal is suppressed by a factor of 100 compared to 
evaporation, practically eliminating the annihilation be¬ 
low ~ 2 GeV. Therefore the spectrum of the evaporating 
particles predicted here does not depend on the nature 
(asymmetric or symmetric) of DM. 

I should also mention that the evaporating DM parti¬ 
cles can create an annually modulated signal with a dif¬ 
ferent phase from the one of the halo DM particles. Here, 
the modulation is due to the small changes in the distance 
of the earth to the sun between summer and winter. The 
perihelion (shortest distance) takes place around January 
3 and the aphelion (largest distance) around July 4. The 
small fluctuation in the distance creates a fluctuation in 
the DM flux arriving on earth, thus the annual mod¬ 
ulation. The largest signal should be expected around 
January 3. Therefore there is a phase difference by al¬ 
most a TT with respect to the annual modulation of the 
halo DM particles. 

Finally, one can study the same effect from evaporating 
DM particles from the earth. Although the capture rate 
of the earth is on average smaller by at least eight orders 
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of magnitude, this can be counterbalanced by the fact 
that the flux is inversely proportional to the distance, 
thus earth is favored by a factor of Additionally, 

the mass below which evaporation dominates is not much 
different from the case of the sun. Although acrit for the 
earth might be a bit smaller compared to the sun, as 
well as the thermalization time slightly larger, the flux of 
evaporating DM particles is not significantly lower than 


the one coming from the sun. Despite this fact, there is a 
fundamental difference. The spectrum of the evaporating 
DM from the earth is dominated by the velocities close to 
the earth’s escape velocity ~ 11 km/sec, which is small 
to create significant recoil for light DM. This is why 1 do 
not examine evaporating DM from the earth here. 
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